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Abstract 

Respiration in bacteria involves a sequence of energetically-coupled electron and proton transfers 
creating an electrochemical gradient of protons (a proton-motive force) across the inner bacterial 
membrane. With a simple kinetic model we analyze a redox loop mechanism of proton-motive force 
generation mediated by a molecular shuttle diffusing inside the membrane. This model, which 
includes six electron-binding and two proton-binding sites, reflects the main features of nitrate 
respiration in E. coli bacteria. We describe the time evolution of the proton translocation process. 
We find that the electron-proton electrostatic coupling on the shuttle plays a significant role in the 
process of energy conversion between electron and proton components. We determine the conditions 
where the redox loop mechanism is able to translocate protons against the transmembrane voltage 
gradient above 200 mV with a thermodynamic efficiency of about 37%, in the physiologically 
important range of temperatures from 250 to 350 K. 

PACS numbers: 87.16.A-, 87.16 Uv, 73.63.-b 



I. INTRODUCTION 



Diffusion-controlled electron and proton transfer reactions are pivotal for the efficient 
energy transformation in respiratory chains of animal cells and bacteria. During the process 
of respiration the energy, extracted from sunlight or from food molecules, is converted into 
an electrochemical gradient of protons (also called a proton-motive force) across an inner 
mitochondrial or bacterial membrane [l, [2], [], [4]. Thereafter, this energy is harnessed by 
ATP synthase for a synthesis of adenosine triphosphate (ATP) molecules, the main energy 
currency of the cell. The energy stored in the proton gradient can be also used to drive a 
rotation of a bacterial flagellar motor. 

The energetically uphill translocation of protons is accomplished by a set of membrane- 
embedded proton pumps or by a redox loop mechanism proposed in the original formulation 
of chemiosmotic theory g For a true proton pump (e.g., cytochrome c oxida^ electro- 
genie events are associated with charges of protons crossing the membrane [2|, |3J. In the 
redox loop mechanism the transmembrane voltage is generated by electron charges moving 
across the membrane. This mechanism is responsible for a proton-motive force generation 
in the respiratory chain of anaerobically grown bacteria such as the facultative anaerobe 
Escherichia coli. In the absence of oxygen and in the presence of nitrate, E. coli can switch 
from oxidative respiration, which uses oxygen molecules as terminal electron acceptors, to 
nitrate respiration, where nitrogen plays the role of a terminal acceptor of electrons in the 
process of nitrate-to-nitrite reduction. 

The redox loop is formed by the formate dehydrogenase- N (Edh-N) enzyme and by the 
nitrate reductase enzyme (Nar) (Fig. 1). The structures of these enzymes and positions of 
all redox centers have recently been determined 6, 7, ^, 9]. As a result of formate reduction, 
HCOO~ — > CO2 + H + + 2e~, a pair of high-energy electrons are delivered to the beginning 
of the pathway (source S) at the P-side of the inner (or plasma) membrane of E. coli. 
Through the intermediate iron-sulfur clusters electrons are transferred, one after another, 
to the integral membrane subunit of Edh-N, which includes hemes bp (site 1) and be (site 
2) located on the opposite sides of the membrane (see Fig. 1). The subindices P and C here 
refer to "Periplasm" and "Cytoplasm", respectively. 

E. coli utilizes a molecule of menaquinone (MQ) as a movable shuttle connecting the 
Edh-N and Nar enzymes. Near the N-side of the membrane menaquinone is populated with 
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two electrons donated by heme be- In this process, menaquinone accepts two protons from 
the N-side of the membrane turning into the form of menaquinol (MQH2). The neutral 
menaquinol molecule diffuses to the P-side where it donates two electrons to heme bi of the 
nitrate reductase and, simultaneously, two protons to the P-side proton reservoir. 

Electrons are transferred, one by one, through heme bi (site 5), to heme bn (site 6) 
and, subsequently, through several iron-sulfur clusters, to the site D on the cytoplasmic 
(N) space where the electrons reduce nitrate to nitrite, NOj — > NOg + H2O. The L and 
H subindices in the notations, b^ and bn, for the sites 5 and 6 refer to "low" and "high" 
redox potentials, respectively. Note that at the beginning of the electron transport chain 
(ETC), where formate is oxidized to CO2 and H + , the midpoint redox potential is very 
low, E m = —420 mV. Thus, electrons entering ETC have high energies (~ 420 meV). 
The menaquinone/menaquinol pair MQ/MQH2 has a much higher redox potential, E m = 
—80 mV (and energy of order + 80 meV), which makes possible the electron translocation 
against the transmembrane voltage. In the second half of the redox loop, formed by nitrate 
reductase, electrons also move energetically downhill, from quinol (E m = —80 mV) to the 
nitrate reduction site having a midpoint potential, E m ~ +420 mV (and energy ~ —420 
meV) 0. 

A geometrical disposition of the quinone-reducing center be and the quinol-oxidizing 
center bi on opposite sites of the membrane is crucial for the generation of the proton- 
motive force |3j, la, LZl] . Electrogenic events resulting in net charge translocation occur when 
an electron moves from heme bp to heme be in the Fdh-N enzyme, and from heme bi to 
heme bn located on the Nar enzyme. 

The crystal structures of the Fdh-N and Nar enzymes solved in Refs. jfj provide key 
components for understanding the mechanism of proton-motive force generation through the 
redox loop. It should be emphasized, however, that the proton-motive force generation is a 
dynamical process, so that structural analysis should be complemented by kinetic studies. 
For example, real time investigations of electron and proton transfers in Complex I 11] and 
Complex IV [12] of mitochondria allow elucidation of a time sequence of transfer events and 
get important information about electron and proton transition rates. Kinetic models of the 
proton pumping processes in cytochrome c oxidase Q] arid in bacteriorhodopsin [ijj] 
are also proven to be beneficial for understanding experimental findings, as well as for an 
initiation of new experiments, giving a comprehensive picture of the phenomenon. 
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In the present work we investigate a redox loop mechanism of a proton-motive force 
generation across the inner membrane of E. coli bacterium within a simple physical model 
incorporating two hemes, bp and be, in the Fdh-N enzyme, two hemes, bi and bjj, in the 
Nar enzyme, and a molecular shuttle (menaquinone) diffusing between these two halves 
of the redox loop. This diffusion is governed by a Langevin equation. There is a pool 
of menaquinone/menaquinol molecules in the bacterial plasma membrane 14, but 

we only consider the contribution of a single menaquinone molecule to the electron and 
proton translocation process. Because of this, the actual values of the electron and proton 
fluxes should be higher than the values calculated below. In order to describe the process 
of loading/unloading the shuttle with electrons and protons, we employ a system of master 
equations, with position-dependent transition rates between the shuttle and electron/proton 
reservoirs. With these equations we analyze the time dependence of the proton-motive force 
generation process together with the dependence of numbers of transferred electrons and 
protons on a transmembrane voltage and on temperature. A thermodynamic efficiency of 
the proton translocation across the inner bacterial membrane is defined and calculated as 
well. 

The paper is organized as follows. In Sec. II we introduce a model of the system and 
present a set of master and Langevin equations, which govern the time evolution of a proton 
translocation process. Sec. Ill is devoted to a discussion of the key parameters of the model. 
In Sec. IV we report our main results and describe the steps for the kinetics of electron and 
proton transfer steps. The conclusions of the paper are presented in Sec. V. 



II. MODEL 



We take into consideration (see Fig. 1) six sites for an electron pathway through the 
system: two sites, 1 and 2, corresponding to hemes bp and be of the Fdh-N enzyme; two 
electron-binding sites, 3 and 4, on the menaquinone shuttle, and two sites, 5 and 6, related 
to hemes bi and bu on nitrate reductase (Nar). For the sake of simplicity we assume that 
heme bp (site 1), located on the periplasmic (P) side of the membrane, is coupled to the 
source of electrons S, and that heme bn (site 6) having a high midpoint potential is coupled 
to the electron drain D. 

The source reservoir S characterized by an electrochemical potential /is and the drain 
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reservoir D described by an electrochemical potential /xd provide a continuous flow of elec- 
trons through the electron transport chain (ETC). The potential fis roughly corresponds 
to the energy of electrons injected into the ETC after formate oxidation, /ig ~ 420 meV, 
whereas the drain potential /iD is related to the electron energy on the nitrate reduction site, 
/io ~ —400 meV. Note that we include the sign of the electron charge in the definition of the 
electron electrochemical potential. This means that a site with a higher electron energy is 
characterized by a more negative redox midpoint potential E m . Here, all energy parameters 
are measured in meV. 

Taking into account two (instead of one) redox sites, 1 and 2, located on opposite sides 
of the membrane, allows us to describe the process of transmembrane voltage generation 
during electron transfer (ET) along the Fdh-N complex. Additional transmembrane voltage 
is generated when an electron moves between two Nar sites, 5 and 6, which are also located 
on the opposite sides of the membrane. 

The pathway for protons includes two proton-binding sites, 7 and 8, on the shuttle. We 
assume that the molecular shuttle moves along a line connecting the redox sites 2 and 5. 
Depending on the position of the shuttle x along this line, the proton-binding sites can be 
coupled either to the positive or to the negative sides of the membrane (P- and N-proton 
reservoirs) . The distributions of protons in the P and N reservoirs are presumably described 
by the Fermi functions with the electrochemical potentials /zp (P-side) and fi^ (N-side of 
the membrane). In its completely reduced form of menaquinol MQH2, the shuttle has a 
maximum load of two electrons and two protons, whereas in its oxidized quinone form 
(denoted by MQ in Fig. 1) the shuttle is empty. 



A. Hamiltonian of the electron-proton system 

Within a formalism of secondary quantization 



13, 
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19 



20J we introduce the creation 
and annihilation Fermi operators, a+, a a , for an electron located on the site a (a = 1, . . . , 6), 
as well as the corresponding Fermi operators, bt,bp, for a proton on the protonable site 
(5 ((3 = 7, 8). The electron population of the a-site is described by the operator n a = a+a Q , 
whereas the proton population of the /3-site has the form: rip = b^bp. Note that we use 
here methods of quantum transport theory to derive classical master equations. A similar 
approach has been applied in studies of quantum coherence in biological systems 21]. 
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The main part of the system Hamiltonian, Hq, involves contributions from the energies, 
e a , of electron sites and energies, e@, of two proton-binding sites on the shuttle complemented 
by terms describing electrostatic repulsions between sites 1 and 2 (with Coulomb energy 
U12) and between sites 5 and 6 (with energy u 56 ). We also add an electron-electron Coulomb 
repulsion between two electron-binding sites, 3 and 4, on the shuttle (with an energy scale 
M34) and a term describing a repulsion between two protons, on the sites 7 and 8, occupying 
the shuttle (energy Uj$). An electrostatic attraction between electrons and protons travelling 
together on the menaquinol shuttle is described by the energy parameters u 37 , u 38 , u 47 , and 
u 48 . As a result, the basic Hamiltonian Hq of the electron-proton system has the form: 

6 8 
a=l /3=7 

+ u 78 n 7 n 8 - u 37 n 3 n 7 - u 38 n 3 n 8 - u 47 n 4 n 7 - u A8 n 4 n 8 

+ (n 3 + n 4 - n 7 - n 8 ) 2 U s (x). (1) 

The last term in Eq. which depends on the shuttle position x, describes the contribution 
of a potential barrier U s (x), which prevents a charged shuttle from crossing the interior of 
the lipid membrane. The barrier has an almost rectangular shape, 

-1 



UJx) = U„ 



x + x 

6X P 1 — j — ) + 1 



(2) 



with a height U s q, a steepness l s , and a width 2x s . This is multiplied by the shuttle charge 
squared: (n 3 + n 4 — n 7 — n 8 ) 2 . The height U s0 of this potential is roughly equal to the energy 
penalty (in meV) for moving a molecule with a charge go (i n units of |e|) and a radius tq (in 
nm) from a medium with a dielectric constant ei to a medium with a constant €2 22j , 



2r \ e 2 ei_ 

For example, the transfer of a charged molecule (go = 1) with radius ro = 0.3 nm, from water 
(ei = 80) to the lipid membrane with e 2 = 3, results in the dielectric penalty U s0 = 770 
meV. The specific shape of the barrier U s (x) in Eq. (j2J) is of little importance for the results 
from this model. 

Electrons in the source (drain) reservoir are described by the creation and annihilation 
operators c^ s ,Cfcs (c^ D ,CkT>), and for protons in the N (P) reservoir we introduce operators 
dq N , dgTss (dgp, d q p), so that the Hamiltonian of the electron source and drain reservoirs, i^sD? 



and the Hamiltonian of the proton reservoirs, H^p, can be expressed as 



k 



1 



(4) 



Here, Eks and Ekr> are the energies of the electrons in the S and D reservoirs, and depend on 
the quasi-momentum parameter k. The energies of the protons in the N- and P-reservoirs, 
£ 9 n and e q p, depend on another continuous parameter q. 

Electrons in the source and drain reservoirs (q = S,D) and protons on the negative and 
positive (a = N,P) sides of the membrane can be characterized by the corresponding Fermi 
distributions, f q (£kq) and F a (e qa ): 



We introduce here the electrochemical potentials n a of the proton reservoirs and the poten- 
tials /x f for the electron source and drain. The potential /is is related to the highest occupied 
energy level of the molecular complex S supplying the ETC with electrons, and the potential 
/id plays a similar role for the molecular complex D providing an electron outflow. 

Couplings between the electron site 1 (heme bp) and the source S, and between the site 
6 (heme &#) and the electron drain D are determined by the Hamiltonian 



with the corresponding transition coefficients tks and tku- The similar Hamiltonian describes 
proton transitions between the shuttle and the proton reservoirs, 



Here, the coefficients T 9 n and T q p, which are assumed to be the same for both sites 7 and 
8, depend on the shuttle position x. The transitions between the redox sites 1, 6 and the 
electron source S and drain D as well as between the N- and P-sides of the membrane and 
the protonable sites 7, 8 on the shuttle are determined by the energy-independent electron 






(6) 




(7) 
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and proton rates mm 



18 



,191,12 



7 S = 27i^2\t k ,\ 2 5{E -e k ,), 

k 

r a = 27iY,\T g *\ 2 S(E-E qa ). 



(8) 



The proton transition rates Tn, Tp depend on the distances (either x + xq or xq — x) between 
the shuttle and the N or P-sides of the membrane: 

1 -2 

Tn = Tno 



x + x 
exp | — ; ) + 1 



PO 



exp 



lp 

Xq — X 

L 



+ 1 



-2 



(9) 



where x = x(t) is the coordinate of the shuttle and l p is the proton transition length. 

The electron tunneling between the redox centers 1, . . . , 6 is governed by the Hamiltonian 

Htmi, 



Htmi — ~ A^a^c^ — A 2 3a^"a3 — A 2 4a 2 f a 4 

- A 35 a^a 5 - A 45 a|a 5 - A 56 a^"a 6 + H.c. 



(10) 



The electrons are transferred between the site 2, located at x = —xq, and the electron- 
binding sites 3 and 4 on the shuttle. On the opposite side of the membrane, at x = xq, 
the electrons tunnel from the sites 3 and 4 to the site 5. These transfers drastically depend 
on the shuttle position x. According to quantum mechanics, we can model the position 
dependence of the tunneling coefficients by the exponential functions: 

\x -\~ Xq 



1 A I 2 - 
1^231 — 


|A 24 | 2 = 


A 2 2 exp | 


-2 


1 A I 2 - 
l^35| — 


|A 45 | 2 = 


|A 5 | 2 exp 


-2 



le 

\x — Xq\ 

L 



(11) 



where l e is an electron tunneling length. 



B. Environment 

The atomic motion of the protein medium has a si gnifi cant effect on electron charge 
transfer between the active sites. Usually (see Refs. 23j, |24j, |25j) the environment is repre- 



sented as a collection of independent harmonic oscillators. The coupling of these oscillators 

8 



to electronic degrees of freedom can be described by the Hamiltonian H env , 

P 2 1 / 6 V 

#env = ^T~. + O m J' W i \ X i X i aUa ~ X i SnS ~~ ^'D^D • (12) 



2m j 2 

3 3 \ « =1 



Here, Xj and pj are the position and momentum of the j-oscillator, having mass rrij and a 
frequency ujj. Also, ns = J2k c ts c ks an d hd = J2k c kD c kr> are the total populations of the 
source and drain reservoirs; Xj a ,Xjs,XjD are the set of coupling constants between electrons 
and their surroundings. 

Thus, the total Hamiltonian of the system has the form 

H = H o + H SD + H NP + H c + Hp + H tun + if e nv (13) 

A unitary transformation, H' = U + HU, with 

U = exp j-iJ^Pj ^2 x jan a + x jS n s + x jD n D j J , (14) 

applied to the Hamiltonian H, removes the environment variables {xj} from the Hamiltonian 
H cnv and introduces phase shifts into the tunneling Hamiltonian H tun : 

H> = H + H SD + H NP + H e + Hp + Hi n + Y,(^ + 1 ^^), (15) 

3 

where 

-^tun = —Qu^l^ — Q23< 1 2 a 3 ~~ QlA&X&A ~ 

Qzh(4, a ^ ~ Q^at a 5 - Q^aia-e + H.c, (16) 
is a new tunneling Hamiltonian, and 

Qaa> = Qa> a = ^aa' exp{j ^ p j ;(x ja - X ja >)}, (17) 

3 

is a phase shift corresponding to the electron transition from site a' to site a (h = 1). For 
simplicity, we neglect here the phase shifts for transitions between the source reservoir and 
the site 1, Xjs = Xj±, and between the site 6 and the electron drain, Xj & = Xj^, together 
with shifts related to proton transfers. The electron and proton reservoirs are described 
by continuous energy spectra. The broadening of the reservoir energy states allows non- 
resonant transitions, e.g., between site 1 and the source S, thus reproducing some effects 
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of the corresponding (1-to-S) phase shifts. Recall also that the tunneling rates A aa i for 
transitions between the sites 2 and 3, 2 and 4, 3 and 5, 4 and 5 depend on the shuttle 
position x(t) and, thus, depend on time (see Eq. (TTTI) ) . However, this time dependence is 
much slower than the time variations of environment-induced phase factors. 

C. Basis states 

To describe all possible occupational configurations of the electron-proton system, we 
introduce a basis of 256 eigenstates, of the Hamiltonian H : H \p) = E^p), /i = 
1,...,256, characterized by the energy spectrum E^. The basis begins with the vacuum 
state, where there are no particles on the sites 1, . . . , 8: |1) = IO1O2O3O4O5O6O7O8), and finally 
ends with the state 1 256) describing the fully-populated system: 1 256) = | I1I2I3I4I5I6I7I8) ■ 
Here, the notation Oq^Iq,) means that the electron site a is empty (occupied). Similar 
notations are introduced for the proton sites 7 and 8. 

It is of interest that all operators of the system, except the operators of the electron and 
proton reservoirs, can be expressed in terms of the basic Heisenberg operators p^ v = 
for example, 

a a = '^^ a a;fMuPfMV, bp — bp.^p^, (18) 

where a, a' — 1, . . . , 6; /3 = 7, 8; and 

a a ;tiv = (fi\a a \u), bp.^ = (n\bp\v) 

are the matrix elements of the electron and proton operators in the basis The Hamilto- 
nian Hq has a diagonal form, 

256 

whereas the tunneling Hamiltonian H tnn (we drop hereafter a prime sign) has only off- 
diagonal elements, 

#tun = - + H.c. (20) 
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Here p M denotes a diagonal operator, p M = p wt = |p)(/i|, and is a combination of 
operators, describing the environment, 



+ Q35{( 1 t a 5)v + < 545(fl4 _ C t 5)^ + Qse^^e, 



fj,v 



(21) 



The Hamiltonian H e , modelling the electron transfer from the source and drain to the sites 
1 and 6, and the Hamiltonian H p , which is responsible for proton transitions between the 
shuttle and the proton reservoirs, are also expressed in terms of the basis matrix p^, 



k fj.v 

H v = -^^(T qN d+ N + T q pd+ p )(b 7 . tlv + b 8 . tlu )p^ + H. 

q fiu 



C. 



(22) 



D. Master equation 



The average value, of the operator p M determines the probability to find the system 
in the state This probability can be found from the Heisenberg equation, 



tun J — j 



(23) 



averaged over the states of reservoirs and over fluctuations of the environment. It is con- 



venient to emr 



systems 



oy methods of quantum transport theory and the theory of open quantum 
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m 



20 



26] to derive the set of master equations describing the time 



evolution of the probability distribution (p M ): 

V V 

where the transition matrix, 



(24) 



V/J, 



k 2a)vii + {^35)ua + {K>4£,)vu + (^56 



IVfll 



(25) 



is represented as a sum o 



the redox states 



24 



27|, 



Marcus rates, (K, aa >) ull , associated with allowed transitions between 



71 



Xaa'T 



(a+a a ') M „| 2 + \(a^a a ') ul j,\ 2 ] exp 



4X aa /T 



,21 



(26) 



11 



where u^v = En, — E Vl and X aa i is the reorganization energy corresponding to the electron 



transition between a to a' redox sites 



18 



20 



The relaxation matrix 7^ describes a 



contribution of transitions between the active sites and the electron and proton reservoirs, 

Ta«/ = 7a{|ai;A«/| 2 [l - /s(^)] + K;^| 2 /s(^v)} + 
7d{| \ 2 [l-hM) + \ I /b(^v)} + 
| 2 )[1 - F N (^)] + (\b 7 , Ufl \ 2 + IVmH^nK,)} + 



(27) 



E. Coulomb energy and redox potential of the shuttle 

The electrostatic coupling between electrons and protons travelling together on the 
menaquinol molecular shuttle is of prime importance for the electron-to-proton energy con- 
version. For the sake of simplicity and without loss of generality, we describe all electrostatic 
interactions involved in Eq. (CQ) by a single electrostatic energy m : M37 = 113$ = U47 = ^48 = 
Uq, and M34 = uj 8 = u . It should be noted that the present model tolerates a significant 
spread (at least 20% and sometimes larger) of the electrostatic parameters. The energy 
scale Un is related to the redox potential E m of the MQ/MQH 2 couple, which is about —80 



meV lOj. To find this relation we model a process of redox titration of a molecule, which 
has one electron and one proton-binding sites characterized by the energy levels e e and e p , 
respectively. 

The electron-binding site is connected to the reservoir of electrons with an electrochem- 
ical potential /x e , whereas the protonable site is coupled to the proton reservoir with an 
electrochemical potential \i v . The energy of the electron-proton Coulomb attraction is de- 
termined by the parameter uq. The goal here is to determine a relation between the electron 
potential \i e and the energy scales e e and u when the electron-binding site is half-populated. 



According to the redox titration procedure 



29] this value of the "ambient" potential (fi 



eJl/2 



determines the redox potential of the molecule E m in the presence of electron-proton elec- 
trostatic coupling, E m = — fi e ,i/2- As in the case of quinone/quinol molecule the protonable 
site should be populated if and only if the electron-binding site is fully occupied. This occurs 
at the condition: 

E p > jJ, p > E p — Uq. 
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Thus, the average electron, (n e ), and proton, (n p ), populations of the molecule are expressed 
in terms of the Fermi distribution function f(e) of the electron reservoir: 

<"> = = i + /(J- 'L-u y < 28 > 

The molecule is half-populated with an electron, (n e ) = 1/2, and with a proton, (n p ) = 1/2, 
when 

^,1/2 = — E m = e e — A (29) 
Calculations for a molecule having two electron sites (with energies e 3 = = e e ) and two 
proton-binding sites (with the energy levels = e$ = e p ) also show the validity of the 
relation Eq. (I29p for the case of a single electrostatic parameter uq. 

F. Proton-motive force 

The difference of proton electrochemical potentials, Afi = fip — fi-^, defines the trans- 
membrane proton-motive force, Afi, consisting of a voltage gradient V and a contribution 
of the concentration difference, ApH, between the sides of the membrane 

A// = V - 2.3 (RT/F) x ApH. (30) 

We introduce here the gas constant R and the Faraday constant F. The potentials 
and V are measured in meV, whereas temperature T is measured in Kelvins (ks = 1). 
At room temperature, T = 298 K, and at the standard gradient of proton concentrations, 
ApH = — 1 , the voltage part of the proton-motive force dominates over the contribution 
of the concentration gradient: Afi ~ V + 60 meV. For example, at Afi — 200 meV the 
voltage difference V ~ 140 meV is applied across the membrane. As a consequence of this, 
the energies, e a , of the redox sites located on the Fdh-N and Nar enzymes are shifted from 
their original values , 



1 



(0) + (_!)- V , (31) 



where (a = 1, 2, 5, 6). We assume here that the voltage drops linearly across the membrane 
[ljj, so that the positions of the energy levels of the electron and proton-binding sites on 
the shuttle are linear functions of the shuttle coordinate x: 

s 7 = e 8 = 4°) + 2^ V, (32) 
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Here, e e and e p are the original values of the electron and proton energies of the shuttle. 



G. Langevin equation 



Within the present model, the Brownian motion of the molecular shuttle 30, |31[ along 
a line, which connects the site 2 (x — — xo) and the site 5 (x — xq), is governed by the 
one-dimensional overdamped Langevin equation 



(x = - 



dUJx) 



((n 3 + n A - n 7 - n 8 ) 2 ) 



dll* (x) 



(33) 



dec doc 
where ( is the drag coefficient of the shuttle in the lipid membrane. The zero-mean val- 
ued, (£) = 0, fluctuation force £ has Gaussian statistics with the correlation function: 
(£,{t)£,{t')) = 2(T5(t — t'), proportional to the temperature T of the environment. The 
diffusion coefficient D of the shuttle is determined by the Einstein relation: D = T/(. The 
potential U c (x), 



U c {x) = U, 



CO 



exp 



I, 



~i -i 



1 



exp 



x + x c 
I, 



1 



(34) 



is responsible for the spatial confinement of the menaquinone/menaquinol molecule inside 
the plasma membrane with the barrier height U c o, the width 2x c (x c > xq) and the steepness 
l c . We also include in Eq. (|33|) the potential U s (x) in Eq. ([2]) hampering the Brownian motion 
of the charged shuttle across the lipid membrane. 



III. PARAMETRIZATION OF THE MODEL 



A. Electron transport chain 

Within our model the electron transport chain begins with the source reservoir S charac- 
terized by the chemical potential /is, which is related (with an opposite sign) to the redox 
energy of formate, fi$ = 420 meV [6|. The redox potentials of hemes bp (site 1) and be 
(site 2) located in formate-dehydrogenase (Fdh-N) are not known. We choose the following 
values: = 445 meV and £3 = 260 meV, for the intrinsic energies of sites 1 and 2. Notice 
that with the transmembrane voltage, V = 140 meV, the energy (see Eq. (13T1) ) of the site 1, 
E\ = 375 meV, is below the potential /is, which is a necessary condition for electron transfer 
from the source reservoir S to the site 1. 
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The original energy of electron-binding sites on the shuttle, ef \ can be related to the 



redox potential E m of the quinone/semiquinone (MQ~/MQ) couple. It is known [32|, [33] 
that the redox energy of the quinone/semiquinone couple is much lower than the potential 
of the quinone/quinol couple. For example, the potential E m for the ubiquinone/ubiquinone 
(UQ/UQH 2 ) couple is about + 60 mV, and the E m for UQ~/UQ couple in aqueous solution 
is of order of —160 mV (4j. For the redox energy of the MQ~/MQ couple, we choose a 
value E m = —215 meV, which is below the known redox energy, E m = —80 meV, of the 
MQ/MQH 2 couple. This means that the energy level of the electron-binding sites is placed 
at = 215 meV. With Eq. (|29l) we obtain a reasonable estimation for the charging energy 
of the shuttle: 

u = 2 (4°) - fi eA/2 ) = 270 meV, 

at yU e ,i/2 = — -E m (MQ/MQH 2 ) = 80 meV. This value of the charging energy u roughly 
corresponds to the electrostatic interaction of two charges located on the opposite sides of 
the menaquinone molecule [34J at a distance ~0.6 nm, provided that the dielectric constant 
e ~ 9. 

We note that at the voltage difference, V = 140 meV, the energy level of the site 2, 
82 = 330 meV, is higher than the level, ef^ + V/2 = 285 meV, of an electron on the 
shuttle located at the N-side. Because of this, electrons can be transferred from site 2 to 
the menaquinone, followed by the proton uptake from the N-side of the membrane. 

The unloading of the fully populated shuttle occurs at the P-side provided the energy 
of the electrons on the shuttle, ef^ — uq — V/2 = —125 meV, exceeds the energy, 5 5 , of 
the site 5. Here, for V = 140 meV, we choose sufficiently low values, £5 = —170 meV and 
Eq = —215 meV, for energy levels of the redox sites 5 and 6 belonging to the second half 
of the redox loop, whereas the original values are e f = -100 meV and 4° } = -285 meV. 
The corresponding redox potentials of these sites differ from the measured redox levels [h]] 
of heme bi '■ E m ~ 20 mV (site 5) and heme bn '■ E m ~ 120 mV (site 6). It is known, 
however, that the redox potentials obtained as a result of equilibrium redox titrations are not 
always applicable for a description of the electron transfer in enzymes, in particular because 
of cooperativity between the redox centers [k3]. This cooperativity can be induced, e.g., by 
electrostatic couplings between the redox sites 1 and 2: w 12 = 20 meV, and between the 
sites 5 and 6: u^q = 20 meV. In the present model the electron transport chain terminates 
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at the drain reservoir characterized (at V = 140 meV) by the energy scale /id = —260 meV, 
which exceeds the energy, — E m = —420 meV, of electrons at the site of nitrate-to-nitrite 
reduction |10l. 



B. Proton pathway 

Protons are loaded on the shuttle at the N-side (x ~ — xq) provided that the shuttle is 
populated at least with one electron. This condition can be met at Sp = uq/2 when the 
energy, Uq/2 — V/2 = 65 meV, of a proton on the shuttle located at x — —Xq, is higher 
than the potential [in, whereas the proton energy level, —u /2 — V/2 = —205 meV, of the 
shuttle, populated with electrons, is below /i^. We take into account electron-electron and 
proton-proton Coulomb repulsions on the shuttle and assume that V = 140 meV, so that 
the total transmembrane proton- motive force, A/z = /ip — /zn, is about 200 meV [35J with 
/xn = —100 meV and /zp = +100 meV. 

Unloading of protons, which occurs at the P-side of the membrane (x ~ xo), is preceded by 
the electron transfer to the site 5. Then, the proton energy goes up, to the level ej 0) + V/2 = 
205 meV, exceeding the potential fip. It should be noted that the present model is robust 
to pronounced variations (As ~ 50 meV) of electron and proton energy levels (see Fig. 3 
later on). 



C. Other parameters 



It is known 



36j that electrons can be transferred between the redox centers in a nanosec- 



ond range. The proton transfer mediated by the hydrogen-bonded chains can occur in 



nanoseconds as well 



23, 



In view of these findings, we choose the following parame- 
ters controlling electron and proton transitions between the reservoirs and the active sites: 
7s — 7d — 0.5/ns, Tn = Tp = 0.05/ns. We assume that all allowed electron transitions 
between the redox sites are determined by the same energy scale A aa i = 8 /xeV. For the 
transition lengths l e and l p involved in Eqs. ( TTTT) we have the values l e = 0.25 nm, 
l p = 0.25 nm. 

The reaction of the environment is described by the set of reorganization energies \ aa i 



lSl2Q|, 



241 ] . which are also assumed to be the same for every pair a, a' : X aa i = A = 100 meV. 
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A similar value of the reorganization energy has been observed in cytochrome c oxidase 39]. 

The Brownian motion of the shuttle is characterized by the diffusion and drag coefficients 
D and C. For the diffusion coefficient we take the value D ~ 3 • 10 -12 m 2 /s, measured 



in Ref. 



40 



41] for ubiquinone (T = 298 K). The drag coefficient Q can be found from 



the Einstein relation, ( = T/D = 1.37 nN-s/m. The potential barrier U s (x) in Eq. 
which impedes the diffusion of the charged shuttle, is characterized by the energy penalty, 
U s o = 770 meV, steepness l s = 0.05 nm, and half-width x s = 1.7 nm. For the potential 
U c (x) in Eq. (l3"4"|) . keeping the shuttle inside the membrane, we choose the height U c0 = 
500 meV, steepness l c = 0.1 nm, and half-width x c = 2.7 nm. The redox sites are located 
at Xq = i 2 nm. On average, the shuttle travels a distance 2xo between sites 2 and 5 in a 
time At = (2xq) 2 / (2D) ~ 2.7 /xs, which is much longer than the time-scales for electron and 
proton transitions to and from the shuttle. 



IV. RESULTS 



To quantitatively describe the kinetics of electron and proton transfers across the mem- 
brane, we numerically solve the system of master equations (T2"4"j) together with the Langevin 
equation (|33|) for a parameter regime, which provides a robust and efficient proton-motive 
force generation, and also roughly corresponds to the menaquinone/menaquinol molecule 
randomly moving inside the bacterial plasma membrane. It should be noted that the present 
model allows significant variations (~20% and sometimes higher) of the parameter values. 

In Fig. 2 we present the time evolution of the electron and proton translocation process 
at T = 298 K, A/x = 200 meV and V = 140 meV. The shuttle starts its motion at x = xq 
(Fig. 2a) and after that diffuses between the membrane borders (shown by two dashed 
red lines at x = ±2 nm). The total electron population, n e = (113) + (714) (continuous 
blue line), and the total proton population, n p = (n 7 ) + (n 8 ) (dashed green line), of the 
shuttle is shown in Fig. 2b. The electron sites 3 and 4 are populated and depopulated in 
concert: (n 3 ) = (rx 4 ) = n e /2. The same relation takes place for the proton sites 7 and 8: 
(717) = (rig) = n p /2. The populations are averaged over the states of electron and proton 
reservoirs as well as over the state of the environment. No averaging over fluctuations of the 
random force in Eq. ( |33l) has been performed in Fig. 2. 

The total number of protons, iVp (dashed green line), transferred by the shuttle from 
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the N- to the P-side of the membrane, and the total number of electrons, Np, (continuous 
blue line) , translocated from the redox site 2 to the site 5 and, finally, to the electron drain 
D, are shown in Fig. 2c. At the beginning of the process (t ~ 0, x ~ — x$) the shuttle 
is rapidly populated with two electrons (n e = 2) and with two protons (n p = 2) taken 
from the N-side of the membrane (/% = — 100 meV). The fully loaded shuttle diffuses and 
eventually reaches (at t ~ 2 /js) the opposite side where the electrons are transferred to 
the redox site 5 (Nn = 2), and two protons (iVp = 2) are translocated energetically uphill, 
to the P-side of the membrane (fip = 100 meV). Accumulation of protons on the positive 
side of the membrane results in a generation of the proton-motive force. The empty and 
neutral quinone molecule diffuses back, to the N-side of the membrane (Fig. 2a), and the 
process starts again. Notice that, as a consequence of the stochastic nature of the process, 
the proton population n p can be a little bit smaller than the electron population n e of the 
shuttle (see Fig. 2b). The resulting tiny charge makes more difficult for the shuttle to cross 
the potential barrier U s (x) in Eq. (j2j). 

It is evident from Fig. 3 and Fig. 4 that the physical mechanism of proton-motive force 
generation described above tolerates significant variations of system parameters such as the 
transmembrane voltage V and temperature T. In Fig. 3 we show the number of protons, 
Np, translocated across the membrane and the number of electrons, Np,, transferred from 
the site 2 to the site 5 as functions of the transmembrane voltage V at T = 298 K. Each 
point in Figs. 3 and 4 is a result of averaging over 10 realizations. Every realization has a 
duration of 100 /is. We calculate the standard deviations for the number Np of transferred 
protons, ap = a/ (Np) — (Np) 2 , and show these deviations as the error bars in Figs. 3 and 
4. The uncertainty o~p, in the number Np, of translocated electrons is close to the value of 
crp. We choose here a symmetric configuration of the proton electrochemical potentials, 



where the potentials /iN,/^p, and the voltage V are measured in meV, and the temperature 
T is measured in Kelvins. 

It follows from Fig. 3 that this redox loop is able to translocate more than 240 protons 
in one millisecond against the transmembrane voltage V ~ 200 meV, which corresponds to 
the proton-motive force Afi ~ 260 meV. In this case (when Np ~ 265, Np, ~ 270, /xp = 
— /iN = 130 meV, /is = 420 meV, and fip, = —260 meV) the thermodynamic efficiency r\ of 




(35) 
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the energetically uphill proton translocation, 

N P up- fi N 

reaches the value t] ~ 37%. 

We note that, despite the dielectric penalty of 770 meV for a charged shuttle, the average 
number of transferred electrons Ad slightly exceeds the number of protons iVp. Interest- 
ingly, both numbers, N P and Ad, have small dips at V = 140 meV. With increasing the 
transmembrane voltage, V > 280 meV, the electron transport from the site 1 (ei = 305) to 
the site 2 (e 2 = 400), and from the site 5 (£5 = —240) to the site 6 (eq — —145, all energies 
in meV) become energetically unfavorable. As a result of this, the numbers of electrons, A^d, 
and protons, N P , translocated across the membrane drop significantly at high voltages. 

The temperature dependence of the average numbers of protons, Np, and electrons, Ad, 
conveyed by the shuttle is presented in Fig. 4 for V = 140 meV. The system demonstrates 
stable performance with Np ~ 220 protons/ms in a window of temperatures from 250 K up to 
350 K. The initial increase of Np and Nd with temperature is probably due to the fact that in 
a warmer environment the shuttle travels more frequently between the sides of the membrane 
transferring more electrons and more protons. Loading (unloading) the shuttle with protons 
follows its loading (unloading) with electrons. At high temperatures menaquinone spends 
less time in the loading zone (at x ~ —x ), and protons have less opportunity to populate 
the shuttle. Therefore, the gap between the numbers of transferred protons and electrons 
widens with increasing temperature. This means that at high temperatures the shuttle has 
more chances to carry a charge, which obstructs the shuttle's diffusion across the membrane. 
Besides that, at sufficiently high temperatures electrons have not enough time to be loaded 
on the shuttle. A combination of these two features results in the high-temperature decline 
of electron and proton flows shown in Fig. 4. 



V. CONCLUSIONS 



Using a simple kinetic model we have examined the process of proton-motive force gen- 
eration across the bacterial plasma membrane. This model is applied to the redox loop 
mechanism of the nitrate respiration in E. coli. This approach includes two redox sites in 
the first half of the redox loop, two redox sites in the second half, and the Brownian shuttle 



19 



diffusing between the negative (N) and positive (P) sides of the membrane. We show that 
the Coulomb attraction between electrons and protons travelling on the shuttle plays an 
essential role in the energetically-uphill proton translocation from the N-side to the P-side 
of the membrane and, thus, in the proton-motive force generation. We have derived and 
numerically solved a set of master equations, which quantitatively describes the process 
of loading and unloading the shuttle with electrons and protons, along with a stochastic 
Langevin equation for the shuttle position. Our model is able to explain the generation of 
the proton-motive force up to 300 meV in the physiologically relevant range of temperatures 
from 250 to 350 K with a peak thermodynamic efficiency of about 37%. A sequence of 
electron and proton transport events and main characteristics of the redox loop mechanism 
calculated in the present paper can be measured in future experiments aimed on a kinetic 
analysis of the nitrate respiration process in bacteria. 
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FIG. 1: (Color online) Schematic diagram of the redox loop. High-energy electrons are delivered 
from the source S to a redox center 1 (heme bp) located near the periplasmic (P) side of the 
membrane. After that, electrons are transferred across the membrane to a redox site 2 (heme 
be) on the cytoplasmic (N) side. At the N-side two electrons reduce a molecule of menaquinone 
MQ, which also takes two protons from the N-side turning into a molecule of menaquinol MQH2. 
The menaquinone shuttle has two electron-binding sites, 3 and 4, and two protonable sites, 7 and 
8. The neutral quinol molecule MQH2 diffuses freely to the P-side of the membrane, where its 
electron cargo is transferred to the redox site 5 (heme and, via the center 6 (heme bn), to the 
drain D on the cytoplasmic side. The oxidation of the quinol molecule MQH2 by the center 5 is 
accompanied by a release of two protons to the P-side of the membrane. Formate dehydrogenase 
(Fdh-N, with centers bp and be) reduces the quinone molecule MQ. Nitrate reductase (Nar, with 
centers bi and bfj) oxidizes the quinol molecule MQH2. Both of these {Fdh-N and Nar) form the 
redox loop, generating a proton-motive force across the membrane. 
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FIG. 2: (Color online) (a) Time dependence of the position x (in nm) (blue continuous curve) 
of the shuttle, diffusing between the walls of the plasma membrane located at x = ±xq (two red 
dashed horizontal lines), where xo = 2 nm; (b) the total proton (n p = {nj) + (rig), blue continuous 
curve) and electron (n e = (713) + (n^), green dashed curve) populations of the shuttle versus time 
(in /us); (c) the number of transferred protons (Np, blue continuous curve) and the number of 
translocated electrons (iVb, green dashed curve) versus time at V = 140 meV, Afj, = 200 meV, and 
at T = 298 K. Notice that the shuttle is loaded near the N-side of the membrane, at x ~ — xq, 
and unloaded at the P-side, at x ~ + xq. It follows from (c) that the process of shuttle unloading 
is accompanied by a stepwise increase of the number of protons, Np, translocated to the P-side of 
the membrane, and the number of electrons, Np,, transferred to the site 5 and, finally, to the drain. 
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FIG. 3: (Color online) The number of protons, Nj> (blue continuous curve), translocated ener- 
getically uphill, from the N-side to the P-side of the membrane, and the number of electrons, Nn 
(green dashed curve), transferred from the site 2 on the Fdh — N enzyme to the site 5 belonging 
to the Nar enzyme, as functions of the transmembrane voltage V at T = 298 K. In Figs. 3 and 
4 the results are averaged over 10 realizations. Each realization has a time span of 100 /is. Error 
bars (standard deviations) are shown for the number Np of translocated protons. 



25 



200 300 400 

Temperature (K) 



FIG. 4: (Color online) Temperature dependence of the numbers of protons, iVp (blue continuous 
curve, with error bars) and electrons, Nd (green dashed curve), transported by the shuttle at the 
transmembrane voltage V = 140 meV. 
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